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Abstract — In Part I of this paper, we presented a mathematical 
model for communication subject to both network interference 
and noise, where the interferers are scattered according to a 
spatial Poisson process, and are operating asynchronously in 
a wireless environment subject to path loss, shadowing, and 
multipath fading. We determined the distribution of the aggregate 
interference and the error performance of the link. In this second 
part, we characterize the capacity of the link subject to both 
network interference and noise. Then, we put forth the concept 
of spectral outage probability (SOP), a new characterization of the 
aggregate radio-frequency emission generated by communicating 
nodes in a wireless network. We present some applications of the 
SOP, namely the establishment of spectral regulations and the 
design of covert military networks. The proposed framework 
captures all the essential physical parameters that affect the 
aggregate network emission, yet is simple enough to provide 
insights that may be of value in the design and deployment of 
wireless networks. 

Index Terms — Spatial distribution, Poisson field, aggregate 
network emission, channel capacity, spectral outage, stable laws. 



I. Introduction 

THE application of the spatial Poisson process to cellular 
networks was investigated in HI, and later advanced in 
El-El- However, these studies focus mostly on error perfor- 
mance metrics, and do not attempt a characterization of the 
channel capacity and interference spectrum. Furthermore, they 
often ignore random propagation effects (e.g. shadowing and 
fading) (T); assume perfect synchronization between different 
interferers at the symbol or slot level Q; or restrict the 
node locations to a disk in the two-dimensional plane |4), 
0, which complicates the analysis and does not provide 
useful insights into the effects of network interference. In 
E), G), the authors analyze coexistence issues in narrowband 
and ultrawideband networks, but consider only a small, fixed 
number of interferers. 

In Part I of this paper 0, we introduced a framework 
where the interferers are scattered according to a spatial 
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Poisson process, and are operating asynchronously in a wire- 
less environment subject to path loss, shadowing, and multi- 
path fading ll9ll- lfT2ll . Under this scenario, we determined the 
statistical distribution of the aggregate interference, and the 
corresponding error performance of the link. In this second 
part, we characterize the capacity of the link subject to both 
network interference and noise. Then, we put forth the concept 
of spectral outage probability (SOP), a new characterization 
of the aggregate radio-frequency (RF) emission generated 
by communicationg nodes in a wireless network. Lastly, we 
quantify these metrics as a function of important system pa- 
rameters, such as the signal-to-noise ratio (SNR), interference- 
to-noise ratio (INR), path loss exponent of the channel, and 
spatial density of the interferers. Our analysis easily accounts 
for all the essential physical parameters that affect the ag- 
gregate network emission. Furthermore, the concept of SOP 
can be used (e.g. in commercial or military applications) to 
evaluate and limit the impact of network interference on any 
given receiver operating in the same frequency band. 

This paper is organized as follows. Section[Tj]briefly reviews 
the system model introduced in Part I. Section [III] analyzes 
the channel capacity of the system, and presents numerical 
examples to illustrate its dependence on important network 
parameters. Section [IV] derives the PSD of the aggregate inter- 
ference, introduces the concept of spectral outage probability, 
and provides numerical examples of both metrics. Section [Vl 
summarizes important findings. 

II. Model Summary 

We briefly review the model introduced in Part I. As shown 
in 1H Fig. 1], we consider the interfering nodes to be spatially 
scattered in the two-dimensional infinite plane, according to 
a homogeneous Poisson process with density A (in nodes per 
unit area). The random distance of interfering node i to the 
origin is denoted by i?;. For analytical purposes, we introduce 
a probe link which is composed of two nodes: the probe 
receiver (located at the origin), and the probe transmitter 
(node i = 0). 

In terms of transmission characteristics, we consider that 
all interfering nodes employ the same two-dimensional mod- 
ulation and transmit at the same power P. For generality, 
however, we allow the probe transmitter to employ an arbi- 
trary two-dimensional modulation and arbitrary power Po, not 
necessarily equal to that used by the interfering nodes. We 
consider that all nodes employ the same symbol rate 1/T, 
but the signal received from node i is shifted by a random 
delay D it where A - U(0, T)Q The probe receiver performs 

'We use U(a, b) to denote a real uniform distribution in the interval [a, b\. 
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coherent demodulation of the desired signal using a conven- 
tional in-phase/quadrature (IQ) detector. 

The wireless propagation channel introduces path loss, log- 
normal shadowing, and multipath fading. Specifically, the 
overall effect of the channel on node i is accounted for 
by the random phase fa ~t/(0,27r), and the amplitude fac- 
tor ^ b — ■ The term accounts for the path loss; a, is due 
to the multipath fading, and has an arbitrary distribution with 
E{af } = 1; and e aGi is due to the log-normal shadowing, 
with d — W(0, 1)@ 

In the rest of the paper, we consider the scenario where the 
location {Ri}°Z 1 and shadowing {Gi]°^ l of the interferers 
(succinctly denoted by V), as well as the shadowing Go 
affecting the probe transmitter, remain approximately constant 
during the interval of interest. This models a quasi-static 
scenario where the movement of the nodes during the interval 
of interest is negligible. In such case, we condition the analysis 
on V in order to derive a capacity outage probability and a 
spectral outage probability, which are more meaningful than 
the corresponding P-averaged metrics0 Other fast-varying 
propagation effects, such as multipath fading due to local 
scattering, are averaged out in the analysis. 

III. Channel Capacity 

In Part I of this paper, we focused on error performance 
metrics. We now build on the results of Part I and analyze 
the capacity of the link between the probe transmitter and 
probe receiver in [H Fig. 1], subject to aggregate network 
interference and additive white Gaussian noise (AWGN). Un- 
like the simple AWGN channel, here the capacity depends 
on the information available about the channel at the probe 
transmitter and receiver. As in Part I, we assume that the probe 
receiver can perfectly estimate the fading (ao and <f>o) affecting 
its own link, thus ensuring that coherent demodulation of the 
desired signal is possible. The probe transmitter, on the other 
hand, is not able to estimate the channel. This corresponds to 
the scenario where the receiver has perfect knowledge of the 
channel side information (CSI). 

A. Capacity Outage Probability 

We start with the complex baseband characterization of 
the probe link, obtained in Part I by projecting all signals 
onto a cosine-sine orthonormal set. Thus, the complex channel 
output Z can be written as 

oCrGo . 



Z = 



-s + w, 



(1) 



where S is the complex channel input, and W is the combined 
aggregate interference and thermal noise, given by 



W 



E 



= <xG, 



X, 



w. 



(2) 



2 We use E{-} and V{-} to denote the expectation and variance operators, 
respectively. In addition, we use M{pi, a ) to denote a real Gaussian distri- 
bution with mean fi and variance a 2 . 

3 We implicitly assume conditioning on V in the rest of the paper, unless 
otherwise indicated. 



with W ~ M(0,N o )B These are essentially the same 
baseband equations as those given in Part I, except that the 
transmitted constellation symbol a§e> Ba has been replaced by 
a generic input symbol S, with an arbitrary distribution /s(s)- 
This emphasizes the fact that to analyze the channel capacity, 
we need to maximize the mutual information over all possible 
input distributions fs(s), and thus cannot restrict S to belong 
to a specific constellation, such as M-PSK or M-QAM. In 
addition, we impose an average energy constraint on the input 
symbol by requiring that E{|S| 2 } < E s . 

Considering that the interfering nodes are coded and op- 
erating close to capacity, then the signal transmitted by each 
interferer is Gaussian, such that X, ~ 7V C (0, 2V X ) HUH The 
resulting aggregate network interference is thus Gaussian when 
conditioned on V, and the distribution of W in (f2]i is given 
b>0 



~ N c {Q,2AVx + N Q ) 



where 



00 g2crGi 



1=1 
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-21, 



(3) 



(4) 



Note that since A in (|4|i depends on V (i.e., {Ri}°Z 1 and 
{Gi}^), it can be seen as a random variable (r.v.) whose 
value is different for each realization of V . It was shown in 
Part I that the r.v. A has a skewed stable distribution IT4l given 
b>0 



A ~ S (a A = j, Pa = 1, 1A = XnC-y* 2 /" 2 ) , 



where b > 1, and C x is defined as 



G, 



r(2-x) cos(ttx/2) ' 



X = 1, 



(5) 



(6) 



with r(x) = Jq 00 t x ~ 1 e~ t dt denoting the gamma function. 

Because of the conditioning on Go and V, equations ([]}-© 
describe a simple Gaussian channel depicted schematically in 
Fig. 03 The capacity of this energy-constrained, fast fading 
channel with receiver CSI can be written as lfl5l 

G = max I(S: Zlan), 

/ S :E|S|2<B S 

where 7(S; Z|ao) is the conditional mutual information be- 
tween S and Z given a^. The optimal input distribution that 

4 We use A/"c(0, a 2 ) to denote a circularly symmetric (CS) complex Gaus- 
sian distribution, where the real and imaginary parts are i.i.d. jV(0, a 2 /2). 

5 Alternatively, we can follow the same approach as in Part I and argue that 
Xi ~ A/c(0, 2Vx) in a scenario where the interferers employ an arbitrary 
two-dimensional modulation (this is the Gaussian approximation introduced 
m GO Eq. (10)]. In such case, Vx is a function of the constellation of the 
interferers, as shown in (8] Eq. (13)]. 

Y 

6 We use X ~ to denote the distribution of X conditional on Y. 

7 We use S(a, /3, 7) to denote a real stable distribution with characteristic 
exponent a € (0,2], skewness £ [—1, 1], and dispersion 7 £ [0, 00). The 
corresponding characteristic function is 



<f>{w) 



cxp [- 7 M Q (1 - j/3siga(w) tan ^)] , a + 1, 
exp [-7M (1 + j^P sign(io) In |iu|)] , a = 1. 



Figure 1. Channel model for capacity analysis. 
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maximizes the mutual information is therefore Af c (0,Es). 
With this input distribution, 



/(S;Z|oo =5o) =log 2 1 



ale 2crG 



°E S 



rf(2AV x +N ) 



in bits per complex symbol, and thus we obtain the capacity 
of the channel as 



C(G ,P)=E ao <Uog 2 1 



v 2 2(tG 



r?(2AV x + N a ) 



Go, A 



(7) 

in bits per complex symbol, where we have explicitly indicated 
the conditioning of C on the random interferer positions and 
shadowing. For a Rayleigh fading channel, af t is exponentially 
distributed with mean 1 and we can further express (0 in terms 
of the exponential integral function Ei(ar) = — ^T~d , t as 



(8) 




C(G ,P) = 

in bits per complex symbol, where 

e 2aG "E s 



V 



r 2b (2AV x +N ) 



(9) 



is the received signal-to-interference-plus-noise ratio (SINR), 
averaged over the fast fading. 

In the proposed quasi-static model, the maximum rate of 
reliable communication for a given realization of Go and V 
is given by I®-©- Such quantity is a function of the random 
user positions and shadowing, and is therefore random. Then, 
with some probability, Go and V are such that the capacity 
is below the transmission rate R, thus making the channel 
unusable for communication at that rate with arbitrarily low 
error probability. The system is said to be in outage, and the 
capacity outage probability is 



PL = F G( „ V {C(G ,V)<R}, 
or, substituting J8j into (fTUb . 



(10) 



(11) 
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Figure 2. Capacity outage probability P£ ut versus the SNR of the probe 
link, for various interferer-to-noise ratios INR (R = 1 bit/complex symbol, 
A = 0.01 m~ 2 , b = 2, r = 1 m, a i3 = 10). 



B. Numerical Results 

Figures [2] and [3] quantify the capacity outage probability and 
illustrate its dependence on the various parameters, such as the 
signal-to-noise ratio SNR = Eg /No, the interference-to-noise 
ratio INR = E/Nq, and spatial density A of the interferers. For 
simplicity, we consider a case study where all interfering nodes 
transmit equiprobable symbols, belonging to a constellation 
that is symmetric with respect to the origin of the IQ-plane 
(e.g., M-PSK and M-QAM). In this particular case, it is 
shown in [8, Eq. (13)] that V x = E/3, and thus © reduces 
to 

e 2fTG °SNR 

"-^(MlNR + l)- (12) 

To evaluate the corresponding P£ ut , we resort to a hybrid 
approach where we employ the analytical result given in 
dTTb-(fT2b. but perform a Monte Carlo simulation of the stable 
r.v. A according to Ifl6l . Nevertheless, we emphasize that the 
expressions derived in this paper completely eliminate the need 
for simulation of the interferers' positions and waveforms in 
the network, in order to obtain the capacity. 

IV. Spectral Characterization of the Aggregate 
Network Emission 

The spectral occupancy and composition of the aggregate 
radio-frequency (RF) emission generated by a network is an 
important consideration in the design of wireless systems. In 
particular, it is often beneficial to know the spectral properties 
of the aggregate RF emission generated by all the spatially 
scattered nodes in the network. This is useful in commercial 
applications, for example, where communication designers 
must ensure that the RF emission of the network does not 
cause interference to other systems operating in overlapping 
frequency bands. To prevent interference, many commercial 
networks operate under restrictions which often take the form 
of spectral masks, imposed by a regulatory agency such as the 
US Federal Communications Commission (FCC). In military 
applications, on the other hand, the goal is ensure that the 
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Figure 3. Capacity outage probability P£ m versus the transmission rate R, for 
various interferer spatial densities A (SNR = INR = 20 dB, b = 2, ro = lm, 
OdB = 10). 



presence of the deployed network is not detected by the enemy. 
If, for example, a sensor network is to be deployed in enemy 
territory, then the characterization of the aggregate network 
emission is essential for the design of a covert system. 

In Part I and in Section [HI] of this paper, we derived the 
error probability and capacity of a link subject to both network 
interference and thermal noise. We now determine the PSD 
of the aggregate interference process Y(t), measured at the 
origin of the two-dimensional plane in J8] Fig. 1]. The spectral 
characteristics of Y(t) can be inferred from the knowledge of 
its PSD. 

A. Power Spectral Density of the Aggregate Network Emission 

The aggregate network emission at the probe receiver can be 
characterized by the complex baseband random process Y(i), 
defined as 



Y(fi) = YY i (t), 



i=l 



(13) 



where Yi(t) is the received process associated with each 
emitting node i. The signal Yi(t) can in turn be expressed 
for all time t as 

Y<(t) 



RS 



fhi{t,T)X4{t-T)<tr, (14) 



where Xj(t) is the complex baseband transmitted signal, and 
hj(f,r) is time-varying complex baseband impulse response 
of the multipath channel associated with node i. The system 
model described by ( Tl4T > is depicted in Fig. [4] It corresponds 
to a generalization of the model introduced in Part I of this 
paper, where we considered a two-dimensional modulation and 
a fiat Rayleigh fading channel. Since now we are interested in 
analyzing the spectral properties of Y(i), we incorporate in the 
model a generic transmitted waveform Xj(t), not necessarily 
associated with a two-dimensional modulation, as well as a 
generic multipath channel hj(t, r), not necessarily associated 
with fiat Rayleigh fading. Also, since in this section we are 
only interested in the aggregate emission of the network, we 
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Figure 4. Channel model for spectral analysis. 

can ignore the existence of the probe link depicted in (8, Fig. 
1]. In what follows we carry the analysis in complex baseband, 
although it can be trivially translated to passband frequencies. 

In the remainder of this paper, we consider that the 
transmitted signal Xj(t) is a wide-sense stationary (WSS) 
process, such that its autocorrelation function has the form 
R^AhM = E{XJ(ti)X,(i 3 )} = i&(At), where At = 
t2 — t\. We define the PSD of the process Xj(i) as 6>x(/) = 
J r At^/{^x(At)}0 Since different nodes operate indepen- 
dently, the processes Xj it) are also independent for different i, 
but the underlying second-order statistics are the same (i.e., the 
autocorrelation function and the PSD of X, (t) do not depend 
on i). As we will show in the case study of Section HV-Cl if 
Xj(£) is a train of pulses with a uniformly distributed ran- 
dom delay (which models the asynchronism between emitting 
nodes), then it is a WSS process. 

In terms of the multipath channel, we consider a wide-sense 
stationary uncorrelated scattering (WSSUS) channel II 171 - 
l2T1l . so that the autocorrelation function of hj(i, r) can be 
expressed as 

fih i (ii,i2,Ti,r 2 ) = E{h*(i 1 ,ri)h i (t 2) r 2 )} 
= P h (At,T 2 )5(T 2 -n), 

for some function Ph(Ai, r). Such channel can be represented 
in the form of a densely-tapped delay line, as a continuum 
of uncorrelated, randomly-scintillating scatterers having WSS 
statistics. The functions h;(t, r) are considered to be indepen- 
dent for different nodes i, but the underlying second-order 
statistics are the same (i.e., the autocorrelation function of 
hi(t, t) does not depend on i). WSSUS channels are an 
important class of practical channels which simultaneously 
exhibit wide-sense stationarity in the time variable t and 
uncorrelated scattering in the delay variable r. They are the 
simplest non-degenerate channels which exhibit both time and 
frequency fading, and also serve as a good model for many 
radio channels. 

We now wish to derive the PSD of the aggregate RF 
emission Y(t) of the network, and with that purpose we 
introduce the following theorem. 

Theorem 4.1 (WSS and WSSUS Channels): Let h(t, r) de- 
note the time-varying complex baseband impulse response of 
a multipath channel, whose autocorrelation function is given 
by Rh(ti,t2,n,T2). Let u(t) denote the complex baseband 
WSS process which is applied as input to the channel, and 
z(t) denote the corresponding output process of the channel. 

8 We use J 7 X _ >SI {'} to denote the Fourier transform operator, where x and y 
represent the independent variables in the original and transformed domains, 
respectively. 
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1) If the channel h(i,r) is WSS, i.e., i? h (ii, *2, n, r 2 ) = 
i?h(At, n, r 2 ), then the output z(i) is WSS and its PSD 
is given b)0 



outage in that frequency band. This leads to the concept of 
spectral outage probability (SOP), which we denote by -P s ut (/) 
and generally define as 



P s (v,n,T 2 )\ v=f i[S u (f)e j2 * f(T1 - T2 

(15) 

where P s {v,t 1 ,t 2 ) = J 7 At^u{Rh{&t,T 1 ,T 2 )}, and 
5 U (/) is the PSD of u(t). 
2) If the channel h(t,r) is WSSUS, i.e., 
Rh(ti,t 2 ,Ti 1 T 2 ) = Ph(At,T 2 )5(T 2 - ri) for some 
function Ph(At,r), then the output z(t) is WSS and 
its PSD is given by 

S.(f)=Vh(j')\v=fis a (f), (16) 

where 2?h(^) = J P s (i>,r)dT is the Doppler power 
spectrum of the channel h(t,r), and P s (v, r) = 
J~At^v{Ph{At, r)} is the scattering function of the 
channel h(t, r). 
Proof: See Appendix [A] for a proof and an intuitive 
interpretation of the theorem. □ 
In the specific context of ( fT4l ). the theorem implies that 
Yi(t) is WSS and thus the aggregate network emission Y(t) 
is also WSS. Furthermore, the PSD of Yj(t) is given by 

S Yi (/) = ^PM/)*S x (/)], d7) 

where 2?h(/) is the Doppler power spectrum of the time- 
varying multipath channel hi(t,r), and 5x(/) is the PSD 
of the transmitted signal Xi(t). Because the processes Yi(t) 
associated with different emitting nodes i are statistically 
independent when conditioned on V, we can write 



dr\ dr 2 



P:M) = WY(f,V)>m(f)}, 



(20) 



S Y (/) = I> Y *Cf)- 



(18) 



Combining ( TPTI i and ( fT8l , we obtain the desired conditional 
PSD of the aggregate network emission Y(t) as 



Sv(f,V) = A[D h (f)*S x (f)}, 



(19) 



where A was defined in (@J. Note that in $1% we explicitly 
indicated the conditioning of 5y on the random node positions 
and shadowing, V . Since Sy(f,'P) depends on V, it can be 
viewed, for a fixed /, as a r.v. whose value is different for 
each realization of Finally, we recall that A, when seen 
as a r.v., has the skewed stable distribution given in ([5). 

B. Spectral Outage Probability 

In the proposed quasi-static scenario, the PSD of the aggre- 
gate network emission 6>y(/, P) is a function of the random 
node positions and shadowing, V . Then, with some probability, 
V is such that the spectrum of the aggregate emission is too 
high in some frequency band of interest, thus causing an 

9 We use * to denote the convolution operation with respect to variable x. 

l0 S-Y(f,V) is in fact a random process whose sample paths evolve in 
frequency instead of time. For each realization V = Vo, we obtain a sample 
path Sy(/, Po) that is a function of /; for a fixed frequency f = fo, 



where S^{f,V) is the random PSD of the aggregate network 
emission Y(t), and m(f) is some spectral mask determining 
the outage (or detection) threshold at the receiver. The SOP is a 
frequency-dependent quantity and, in the case of slow-varying 
positions V, is a more insightful metric than the PSD averaged 
over V. Note that this definition is applicable in general to any 
emission model: the spectral outage probability f s ut (/) rep- 
resents the probability that the PSD of the aggregate network 
emission, measured at an arbitrary location in the plane and at 
a particular frequency /, exceeds some predetermined mask. 

In commercial applications, the concept of SOP can provide 
a radically different way to establish spectral regulations. 
Current regulations and standards (e.g., FCC Part 15 or IEEE 
802.11) impose a spectral mask on the PSD at the transmitter, 
and the type of mask often depends on the environment in 
which the devices are operated (e.g., indoor or outdoor). The 
purpose of this mask is to limit RE emissions generated 
by a terminal, and to protect other services that operate 
in dedicated bands (e.g., Global Positioning System, public 
safety, and cellular systems). However, the transmitted PSD is 
usually not representative of the aggregate PSD at the victim 
receiver, due to the random propagation effects (multipath 
fading and shadowing) and the random position of the emitting 
nodes. Thus, spectral regulations that are based only on the 
transmitted PSD do not necessarily protect a victim receiver 
against interference. 

The approach proposed here is radically different, in the 
sense that the spectral mask is defined at the victim receiver, 
not at the transmitter. In effect, the mask m(f) introduced 
in (120b represents the outage threshold with respect to the 
accumulated PSD at the receiver, not the individual PSD at 
the transmitter (this follows from the fact that 5y(/, V) is 
measured at an arbitrary location in the plane, where a probe 
receiver could be located). Therefore, the received aggregate 
spectrum Sy(f, V) and the corresponding -P D s ut (/) can be used 
to characterize and control the network's RE emissions more 
effectively, since they not only consider the aggregate effect 
of all emitting nodes at an arbitrary receiver location, but also 
incorporate the random propagation effects and random node 
positions. Furthermore, the use of different masks for indoor 
or outdoor environments is no longer necessary, since the en- 
vironment is already accounted for in our model by parameters 
such as the amplitude loss exponent b, the spatial density A 
of the emitting nodes, and the shadowing coefficient a. 

In military applications, on the other hand, the goal is 
to ensure that the presence of the deployed network is not 
detected by the enemy. If, for example, a surveillance network 
is to be deployed in enemy territory, then the characterization 
of its aggregate emission is essential for the design of a 
covert network with low probability of detection. In such 
application, the function m(/) in d20b can be interpreted as 
the frequency-dependent mask which determines the detection 
threshold (not the outage threshold as before). In other words, 
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if the aggregate spectral density Sy(f, V) measured at a given 
location exceeds the mask m(f), then the presence of the 
deployed network could be detected by the enemy. 

For the signal model considered in this paper, P^ ut (f) can 
be derived by substituting ( fT9b into the general definition of 
SOP in (|20]i, leading to 



PL(f) 



A > 



I -F A 



m(f) 



2M/)*«Sx(/) 
Ph(/)*5 X (/) 



(21) 



where Fa{-) is the cumulative distribution function (c.d.f.) of 
the stable r.v. A, whose distribution is given in (O. 



C. Numerical Results 

We now present a case study to quantify the spectral den- 
sities and outage probabilities derived in the previous section. 
We also illustrate their dependence on the various parameters 
involved, such as the transmitted pulse shape, spectral mask, 
transmitted power, and spatial density of the emitting nodes. 
For all numerical examples, we consider that the emitting 
nodes employ a two-dimensional modulation (e.g., A/-PSK or 
A/-QAM), such that transmitted signal X^(t) can be written 
for all t as 



Xi(t) = ai,n9{t-nT-Di), 



(22) 



where the sequence {&i,n}t=-oo represents the stream of 
complex symbols transmitted by node i, assumed to be in- 
dependent identically distributed (i.i.d.) in n and zero-mean, 
for simplicity; g(t) is a real, baseband, unit-energy shaping 
pulse, defined for all values of t; T is the symbol period; 
and Di ~W(0, T) is a random delay representing the asyn- 
chronism between different emitting nodes. The type of con- 
stellation employed by the emitting nodes is captured by the 
statistics of the symbols {a, i, n }o Note that the process X,(i) 
in (23 is WSS, as required by Theorem |4TF| T 
Xi (i) is then given by 



The PSD of 



5x(/) = J P|G(/)| 2 , 



(23) 



where P = E{|a i n | 2 }/T is the power transmitted by each 
emitting node, and G(f) = J r {.g(i)}. 

In terms of the multipath channel, we consider for simplicity 
that h(i, t) is time-invariant such that it does not introduce 
any Doppler shifts, i.e., T>h(v) = 5(^)0 Substituting the 

"Note that each complex symbol aj n = ai jn e :,8i '" can be represented 
in the IQ plane by a constellation point with amplitude dj n and phase 

12 This can be shown in the following way: first, if we deterministically 
set Di to zero in I22t . the resulting process X;(i) is wide-sense cyclosta- 
tionary (WSCS) with period T (22); then, since Xj(t) = Xj(t- Di), with 
Di ~ W(0, T) and independent of everything else, it follows that X;(t) is 
WSS. 

13 For typical node speeds or channel fluctuations, the frequencies of the 
Doppler shifts are on the order of few KHz. As a consequence, when the con- 
sidered Xj(t) is an ultrawideband signal, Ti^ly) can be well approximated 
by a Dirac-delta function. 
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(a) PSD of the individual transmitted signal versus frequency (bottom curves), 
for various pulse shapes g(t). The square and Hanning pulses are normalized 
so that the transmitted signals have the same power P. The piecewise-constant 
spectral mask m(f) (top curve) determines the outage threshold at the receiver. 
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(b) Spectral outage probability P^utif) versus frequency, for the piecewise- 
constant mask m(f) shown in (a). 

Figure 5. Effect of the transmitted baseband pulse shape g(t) on the PSD and 
the outage probability P s ul (/) (P = lOdBm, T = 1CT 6 s, A = 0.1m- 2 , 
b = 2, a dB = 10). 



expressions for <Sx(/) and T>^{v) in (l2T1 l. we obtain the SOP 

as 

m(f) 



PL(f) = 1-F A 



P\G(fW 



(24) 



Figure [5] shows that for a fixed spectral mask m(f), 
the SOP can be highly dependent on the pulse shape g{t) 
(e.g., square or Hanning pulse). In fact, P% ut (f) is a nonlinear 
function of |G(/)|, where the nonlinearity is determined in 
part by the c.d.f. Fa(-) of the stable r.v. A, as shown in 
d24l . Thus, the SOP can be used as a criterion for designing 
the pulse shape: for example, we may wish to determine 
the baseband pulse g(t) and transmitted power P such that 
max/ Pgatif) < P* i where p* is some target outage probability 
which must be satisfied at all frequencies. 

Figure [6] shows that for a fixed pulse shape g(t), 
P D s ut (/) can significantly depend on the spectral mask m(f) 
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(a) Plot of various spectral masks m(f) which define the outage threshold at 
the receiver (top curves). Also shown is the PSD of the individual transmitted 
signal versus frequency (bottom curve). 
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(b) Spectral outage probability P \ lt (/) versus frequency, for the various 
masks m(f) shown in (a). 

Figure 6. Effect of the spectral mask shape m(f) on the outage probabil- 
ity P s ut (/) (square g(t), P = 10 dBm, T = 10 -6 s, A = 0.1 m" 2 , b = 2, 
CTdB = 10). 
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Figure 7. Spectral outage probability Po Ut (f) versus frequency, for various 
transmitted powers P (square g(t), T = 10 — 6 s, A = 0.1 m~ 2 , 6 = 2, og^ = 
10, m(f) = -60dBm/Hz). 
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Figure 8. Spectral outage probability P s ut (/), evaluated at / = 0, versus 
transmitted power P, for various spatial densities A of the emitting nodes 
(square g(t), T = 10 -6 s, b = 2, <r dB = 10, m(/) = -60dBm/Hz). 



(e.g., piecewise-linear, Gaussian, or constant mask). Since 
Pout(f) accounts for both G(f) and m(f), it quantifies the 
compatibility of the transmitted pulse shape with the spectral 
restrictions imposed through m(/). 

Figures [7] and [8] illustrate, respectively, the dependence of 
the outage probability P£ m (f) on the transmitted power P and 
spatial density A of the emitting nodes. Specifically, as P or 
A increase, the aggregate network emission becomes stronger, 
and thus P s llt (/) deteriorates at all frequencies, approaching 
the maximum value of 1. 

D. Generalizations 

We now extend the results to an heterogeneous scenario 
with K different networks, where a given emitting node 
belongs to the network k £ {1...K} with probability p^, 
independently of everything else. Using the splitting property 
of Poisson processes l26l . we know the emitting nodes from 



each network k form a spatial Poisson process, which is 
independent of the processes of other networks and has spatial 
density A^ = Xpk- Therefore, we can write the aggregate 
emission from all nodes in all networks as 

K 

Y(t) = J2 Y k(t), (25) 

fc=i 

where Yfc(t) = Y^Li Y k,i(t) is the aggregate emission from 
the individual network k, and 

Y M (*) = / h M (*, tJXjm (t-T)dr, k = l...K, 

H k,i J 

(26) 

where Xfc »(t) and hk,i(t,T) are, respectively, the transmitted 
signal and the impulse response of the multipath channel 
associated with node i from network k. We consider that 
Xfc i(t) and hfc r) are independent in both k and i. Then, 
the aggregate emission (t) is also independent for different 



networks k when conditioned on the positions V, and thus 
<Sy(/) = SfcLi "^y,, (/)• We can generalize ( fl9l ) and write 
the conditional PSD of the aggregate emission Y(t) in this 
heterogeneous scenario as 



K 

E 

fc=i 



A fc [©h fe (/)*5 Xji (/)], 



(27) 



where 2?h fc (/) and 5x fc (/) are, respectively, the Doppler 
power spectrum and the PSD of the transmitted signal as- 
sociated with network k; and the r.v.'s {Ak} are i.i.d. in k and 
given by 

00 a 2aG k i 
i=l 



Rf 



k . i 



In (|27| |, we have explicitly indicated the conditioning of 6>y 
on the random node positions and shadowing V. Since Ak 
depends on V, it can be seen as a r.v. whose distribution is 
skewed stable, 



ola = Pa 



1 



1A = Afc7rC 1/h e 



1 2<T 2 /b' 



(28) 



V. Summary 

This two-part paper investigates a mathematical model 
for communication subject to both network interference and 
AWGN, where the spatial distribution of the nodes is captured 
by a Poisson field in the two-dimensional plane. We specifi- 
cally address the cases of slow and fast-varying node positions, 
as well as homogeneous and heterogeneous networks, in a 
realistic wireless environment subject to path loss, multipath 
fading and shadowing. In Part I, we determined the statistical 
distribution of the aggregate interference at the output of 
a conventional linear receiver, which leads directly to the 
characterization of the error performance (in terms of outage 
and average probabilities). 

In this second part, we characterize the capacity of the 
link when subject to both network interference and noise, 
and derive the PSD of the aggregate RF emission of the 
network. Then, we put forth the concept of spectral outage 
probability (SOP), and describe some possible applications, 
namely the establishment of spectral regulations and the design 
of covert military networks. In particular, the SOP can be 
used as a criterion for designing pulse shapes or controlling 
interference in wireless networks, and as a measure of the 
network's covertness. Our framework clearly shows how the 
aggregate network emission can be characterized in terms of 
important network parameters, thereby providing insights that 
may be of value to the network designer. 

Appendix A 
Derivation of TheoremI4.1I 



The derivation of Theorem 14.11 relies on the general the- 
ory of linear time-varying systems and Bello system func- 
tions 1T7| - ||20| . Let h(f, r) denote a time-varying complex 
baseband impulse response of a multipath channel. When the 



complex baseband process u(t) is applied as input to the 
channel, the output process z(t) is given by the integral 



z(i) = J h(t,r)u(f - T)dr. 



We define the autocorrelation function of the input u(t) 
as R u (ti,t 2 ) = E{u*(ti)u(f2)}, and the autocorrelation 
function of the channel h(t, r) as Rh(ti,t2,Ti,T2) = 
E{h*(ii,Ti)h(t2,T2)}. The autocorrelation of the output z(i) 
is generally given by 



R z (tl,t 2 ) = / Rh(tit2,Ti,T2)R u (t 1 -Tx,t2-T 2 )dT 1 dT2. 

(29) 

Since the input process u(t) is WSS, R u (ti,t 2 ) = Ru(At), 
where At = t 2 — t-y. 

We first consider a WSS channel h(t, r) such that 
Rh(tii t2, ti, t 2 ) = Rh(A.t, T\ 1 T2). Then, we can rewrite 
R z (h,t 2 ) in (123 as 



R z (h,t 2 ) = j I R h (At, n ,T2)R u (At + n-T 2 )dT 1 dT2 
= R,(At). (30) 

Since R z (ti,t2) is a function only of At, the output z(i) is 
also WSS. The PSD of z(t) can be written as 

S z (.f) = T At ^ f {R z (At)} 

J R h (At, Tl ,T 2 )R u (At + n - t 2 ) 



x e 



d(At) 



dr\dT2 



P a (v,n,T2)\v=f* Su{f)e 



j27r/(-ri-T 2 ) 



dTidr 2 , 
(31) 



where P s {v, n,T 2 ) = T A t^{Rh(At,Ti,T 2 )}, and S u (f) is 
the PSD of u(t). This is the result in Theorem |4~T1 eq. < fT3T >. 

We now further constrain the channel h(t, r) to be WSSUS 
such that Rh(ti,t 2 ,Ti,r 2 ) = Ph(At,T 2 )S(T 2 — ti), for some 
function Ph(At,r). Then, R z (ti,t 2 ) in (l30l can be further 
simplified as follows: 



R z (h,t 2 ) = J J P b (At,T 2 )5(T2-n)R u (At + n-T 2 )dndT2 
= J P h {At,T)R u (At)d T 

= R u (At) J P h (At,r)dT ±R.(At). (32) 
The output z(t) is therefore WSS, and its PSD can be written 



as 



S z (f) = T At ^ f {R z (At)} 

= T A t^ f {Ru{At)} * T At ^ S j J P h (At, r)dr| 

= S u (f)i I ' P s (v,T)\„ =f dT, (33) 

where P b (v,t) = J r At ^ tU {P\ x (At 1 r)} is known as the scat- 
tering function of the channel h(t, t). It provides a measure 
of the average power output of the channel as a function 
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of the delay r and the Doppler shift v. Furthermore, if 
we define the Doppler power spectrum of the channel as 
£>h(» = / Ps{v,T)dT, then CH can be succinctly written 

as 

5 z (/)=5 u (/)*O h (z/)| !/=/ , (34) 

which is the result in Theorem 14. II eq. ( TTol l. 

From d34l , we conclude that S z (f) depends on the Doppler 
power spectrum of the channel, J P s {y,r)dr, but not on its 
power delay profile J P s (v : T)dv. This is intuitively satisfying 
since all delayed replicas of the WSS process u(t) have 
the same PSD. Furthermore, if the channel h(t, r) is time- 
invariant, then T>h(i/) = 8{v) and thus S z (f) = S u (f), 
i.e., the channel does not affect the PSD of the input. On the 
other hand, if the channel is time-varying in such a way that it 
introduces a Doppler shift of /o Hz, then 2?h(f ) = 8(v — fo) 
and thus S z (f) = S u (f — fo), i.e., the output PSD is simply 
the input PSD shifted by fo Hz, as expected. 
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